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Motivated by the problem of weak collective pinning of vortex lattices in high- 
temperature superconductors, we study the model system of a four-dimensional elastic 
manifold with N transverse degrees of freedom (4+7V-model) in a quenched disorder en- 
vironment. We assume the disorder to be weak and short-range correlated, and neglect 
thermal effects. Using a real-space functional renormalization group (FRG) approach, 
we derive a RG equation for the pinning-energy correlator up to two-loop correction. 
The solution of this equation allows us to calculate the size R c of collectively pinned 
elastic domains as well as the critical force F c , i.e., the smallest external force needed 
to drive these domains. We find R c oc 5 P 2 exp(ai/5 p ) and F c oc 8~ 2a2 exp(— 2a.i/5 p ), 
where S p <C 1 parametrizes the disorder strength, a± — (2/ty) n ^ 2 8tt 2 /(N + 8), and 
Q.2 = 2(5iV + 22) /(N + 8) 2 . In contrast to lowest-order perturbation calculations which 
we briefly review, we thus arrive at determining both qi (one-loop) and ot2 (two-loop). 



I. INTRODUCTION 

The generic problem of a D-dimensional elastic manifold with TV transverse degrees of freedom subject 
to quenched random impurities that can pin the manifold has attracted a great deal of attention for many 
years. It has applications in numerous areas of physics such as dislocations, spin systems, polymers, 
charge density waves, and vortex lattices in type-II superconductors. The case N = 1, for instance, 
corresponds to D-dimensional interfaces separating two coexisting phases in (D + l)-dimensional systems. 
On intermediate distances the flux lattice in superconductors is believed to behave like an elastic manifold 
with D = 3 and N = 2. 

The crucial quantity for defining order in these systems is the disorder-averaged square of the relative 
displacement (u 2 (x)) := ([u(x) — u(0)] 2 ), which describes the roughening of the TV-component displace- 
ment field u(x) of £>-dimensional support. In general, (it 2 (x)) depends on the distance R := |x|, the 
dimensionalities D and N, and on the type of disorder. In the present paper, we choose the disorder to 
be weak, short-range correlated, and of Gaussian type. Under these assumptions it has been shown |]| 
that in less than four dimensions (D < 4) quenched disorder always destroys the translational long-range 
order of the manifold, breaking it up into correlated domains of extent R c . Each of these domains behaves 
elastically independently and is pinned individually. The length scale R c over which short-range order 
subsists is commonly referred to as the collective pinning radius and is usually defined by the criterion 



(u 2 (R c )} ~ £ 2 , where £ is the length scale describing the internal structure of the manifold J_j9[]. For 
distances R smaller than R c (perturbative regime), the ground state of the system is unique and the mean 
square displacement (u 2 (R)) can be calculated perturbatively. On the other hand, for R > R c (random 
manifold regime), there are many competing metastable minima inducing the so-called "wandering" of 
the manifold as described by the scaling law (u 2 (R)} « £ 2 (i?/i? c ) 2< > (R ^> R c ). A lot of effort has gone 
into the determination of the wandering exponent £, using elaborate techniques such as the replica for- 
malism combined with either variational approaches §|| or with the renormalization group (RG). 
Here, we make use of RG methods within the perturbative regime to determine the collective pinning 
radius R c . For that purpose, we will replace the estimate (u 2 (R c )) ~ £ 2 by an improved definition of R c 
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based on the divergence of the fourth derivative of the pinning-energy correlator, signaling the appearance 
of competing ground states at this length scale. The collective pinning radius R c is not only a relevant 
quantity for characterizing order, but also determines dynamic properties such as the minimum force F c 
needed to move the pinned manifold (critical force) and the activation barrier U c for creep. 

In this work we concentrate on the particularly interesting case D = 4, the upper critical dimension 
of this problem, where the mean-square displacement shows a logarithmic behavior in the perturbative 
regime JTo|] , (u 2 (i?)) oc lni?, (R < R c ). This logarithmic situation motivates the use of renormalization 
group methods. The (4+7V)-dimensional model system is strongly related to the problem of weak collec- 
tive pinning of a (3+2)-dimensional manifold with dispersive elastic moduli, an issue which has received 
a lot of interest in connection with pinning and creep behavior of vortices in high-T c superconductors. In 
a future paper pO] , the RG technique developed below will be applied to that topic (see also Ref. [19]). 

The main objects of interest in our investigations are the collective pinning radius R c and the critical 
force density F c . In section II, we introduce the non-dispersive (4+A)-model and define the static and 
dynamic Green's functions as well as the relevant correlation functions. For simplicity, we neglect thermal 
effect, i.e., we set T = 0. Then, in section III, we derive R c and F c by means of lowest-order perturbation 
calculations combined with scaling techniques and the dynamic approach, and discuss the problem with 
these methods. Section IV is devoted to the renormalization group treatment of the pinning problem. We 
construct the RG transformation and derive a functional RG equation for the pinning energy correlator 
in one and two-loop approximation. We show how R c can be obtained by solving this equation and 
determine the critical force density F c using simple scaling relations. We discuss various aspects making 
this technique superior to those presented in section III. In a second step, we apply these results to the 
dynamical situation, where the system is driven with an external force. In particular, we determine the 
behavior of the friction coefficient r\ under the action of the RG and give an alternative derivation of F c . 
Finally, in section VI we summarize and present our conclusion. 

II. MODEL 

A. Non-dispersive (4+iV)-model 

We consider a four-dimensional clastic manifold in a (4+A)-dimensional space in the continuum limit: 
A point of the manifold is represented by a four-component vector x (internal degrees of freedom), while 
the displacement relative to the equilibrium position at that point is characterized by the continuous 
TV-component vector field u(x) (transverse degrees of freedom). The elastic part of the free energy is 
given by the isotropic expression 

F el [u] = - C J d 4 xVu a (x) • V«"(x) , 

where a = 1, 2, . . . , N; as usual, indices appearing twice are implicitly summed over. In the present work, 
the elastic modulus C is taken to be a constant (non-dispersive model). 

The elastic manifold is embedded in a random environment which we model as a Gaussian random 
pinning potential C/ p i„(x, s) with zero mean and short-range correlations, 

(<7 pm (x,s)[/ pm (x',s')} = 7 t7 ,5 4 (x-x')<5 JV (s-s'). (1) 

The parameter 7(7 is a measure of the disorder strength and is assumed to be small. Angular brackets 
(. . .) denote the average over all possible realizations of disorder. Note that U P i n depends on the internal 
as well as on the transverse degrees of freedom of the system. For a given configuration u(x) and a 
random, but fixed disorder, the energy density arising from the interaction of the manifold with the 
disordered potential is the convolution of the potential U P i n with the form factor p(s): 

E pin (x,u(x)) = J d N sU pm (x,s)p(\s - u(x)|) . (2) 

The function p(s) characterizes the internal structure of the manifold. In the following, we will adopt the 
simple expression 
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p(s) = exp ^ j , (3) 

where £ is chosen as the smallest resolvable length scale of the system. The effect of the form factor is 
to smear the random potential U P i n over a length £, so that the typical effective distance between two 
'valleys' in the energy landscape is of order £. 

The total free energy of the manifold in the presence of disorder is eventually 



^GW(x) ■ Vu Q (x) + E pin (x, u(x)) 



(4) 



with the pinning energy correlator 

(£ pin (x,u(x))£ piB (x>'(x'))) = <5 4 (x-x')^(u(x)-u'(x')), (5) 

where 



K s (u)= 7u e 2 exp (-^j , u:=\u\, (6) 

as obtained from Eqs. The exponential in the correlator K^(u) is directly related to the expression 

for the form factor, Eq. ([|). 

The equilibrium configuration u(x) of the system results from the competition between the elastic 
energy the system has to pay for a distortion and the pinning energy the system can win by accommo- 
dating to the impurity potential. Unfortunately, determining these configurations reveals an impossible 
undertaking because of the randomness of E pin . One must therefore content oneself with calculating 
disorder-averaged quantities, relying on the system's assumed property of self-averaging. 



B. Green's functions 

Minimizing the free-energy functional ^ with respect to u(x) 
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defines the equation of state, 



^[u] = 0, a=l,...,N, (7) 



where F p ^ n stands for the pinning force density, 



CVV(x) = F; ra (x,u(x)), (8) 



S d 

^pm(x,u(x)) := - ^7T ^[U] = -— £pi n (x,u(x)) . (9) 
The equation of state (|^) can equivalently be written in the form 

u°-(x) = J AG^(x-y)^„(y,u(y)), (10) 



where 



G^(x)=^G(x), G(x) = {2n) l Cx2 , x:=\x\, (11) 



is the static Green's function of the system defined by the relation 

-GV 2 G(x) =5 4 (x). 
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For the dynamic approach in section III we need the time-dependent Green's function G a ^(x, t). We 
assume a dissipative dynamics for our elastic manifold which we characterize by the friction coefficient 
T] . With the equation of motion 

-CV 2 + Vo U Q (x, t) = F«„(x, u(x, t)) , 

the dynamic Green's function takes the form G a @(k, u) — S a ^/(Ck 2 — irj Lo) in Fourier representation. 
Transforming back to real space, one readily obtains 



G * 3 (x, t) = &<**> G(x, t) , G(x, t) = exp (- 



r) x 2 



(12) 



4G£ 

The Heaviside-function O(i) reflects the causality of the Green's function. 

C. Correlation functions 

For later use we define the derivatives of E P i n , 

(x, u(x)) := — • • ■ g^E pin (x, u(x)) , n, n, 1 V. Z = 1, 2, 3, . . . , 

and the corresponding correlators, 

(£; p 7„- Qi (x,u(x))^-^(x',u'(x'))> = (-iri 4 (x-x')^- a ' ft -^(u(x)-u'(x')), (13) 

where naturally (u) :— (d/du a )(d/du 13 ) ■ ■ ■ i<Q(u). Each partial derivative d/du' a contributes a 

minus sign, thus giving rise to the factor (—1)™. In particular, the force-force correlator is given by 

«u) . - 7 „ {»- (=) 4 («- - ^) exp (- * ) . (14) 

Of special interest in section IV will be the curvature of the force at the origin, 

Kf lS {0) = T e (8 af} 5 lS + S^S 135 + S^S 131 ) =: T 6 , (15) 

JV 

with r 4 :=iu^ N - 4 (f) 2 . 

III. PERTURBATION THEORY 

In this section we present the usual ways to define and determine the collective pinning radius R c and 
the critical force density F c . First, we calculate R c by a lowest-order perturbation calculation, and assess 
F c by dimensional estimates |fl]]. Second, we use the dynamic approach |l2|,[l3| to determine F c in an 
alternative manner. We also point out the disadvantages of these simple methods. 

A. Dimensional estimates 

The collective pinning radius R c is usually derived by computing the fluctuations of the displacement 
field u(x) induced by the random environment, 

( M 2 (^):=([u(x)-u(0)] 2 ), R:=\x\, 

in the absence of an external force, F ext = 0. The condition (u 2 (R c )) ~ £ 2 defines the collective pinning 
radius R c which represents the extension of a collectively pinned elastic domain. For small fluctuations, 
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(u 2 (R)} < £ 2 , one can use a perturbative approach jflj. The starting point is the integral equation (10). 
We expand the pinning force, F pin (y, u(y)) = F pin (y,0) + 0(u(y)), and only retain the lowest-order 
term, so that Eq. @ reads (G Q/3 = 6 al3 G) 

u a (x) = J d 4 yG( X -y)F£ n (y,0), (16) 

From this expression and by means of Eqs. (|9|) and (13), we obtain 

(u«(xK(x')) = -K^(0) J d 4 y G(x - y) G(x' - y) + 0(<#) . 

Next, we write the Green's functions in Fourier representation; the mean square of the relative displace- 
ment then takes the form 

(u 2 (R)) = -2^n°) / (! - cos ( k ■ x )) G'(k) 2 + 0(7&) • (17) 

Surprisingly enough, this lowest-order perturbation expression is true to all orders in ju flbH (this can 
be checked by keeping u(y) on the rhs. of Eq. (16) and going through the subsequent calculations). Two 
elements are responsible for this striking fact: the short-range correlation oc <5 4 (x — y) in (E p i n E p i n ) and 
the implicit assumption made in this derivation that the system is in a unique state, which is only valid 
in the perturbative regime. The integral in Eq. (|l^) is invariant under space rotations; one can therefore 
choose x parallel to the a^-axis and introduce cylindrical coordinates (K 2 :— k\ + fc| + k 2 ) , 



Performing the iC-integration yields 

M*)> = -^^/ o dk, ^ • (18) 

The ^-integration produces a log-contribution which diverges for &4 — > cx>. For this reason, one has to 
introduce an upper cutoff given by the inverse of the smallest length scale in the problem which is £. 
With Kg a (0) = -7,7 £ N ~ 2 N(tt/2) n / 2 , Eq. (0, one finally obtains 

where 5 p := 7(7 £ "~^/CP is the dimensionless disorder parameter. The collective pinning radius defined 
through the relation (u 2 (R c )) ~ £ 2 then is 

which is exponentially large in the limit of weak pinning, S p -C 1. With the upper cutoff in Eq. ( |l8|) being 
given only up to a factor of order unity, the constant of proportionality in Eq. (20) is not unequivocally 
determined. Much more relevant, though, is the uncertainty in the numerical factor in the exponential 
function, which has its origin in the criterion (u 2 (R c )) ~ £ 2 , where the constant of proportionality is 
assumed to be of order unity. 

A rough estimate for the critical force density F c can be gained by scaling. The typical elastic energy 
of a collectively pinned domain is 



u °~ c {±) v ~ 

where V c :— R A is the collective pinning volume. The energy gain arising from the action of an external 
force F ext shifting the domain by £ is 

Uext ~ Fext Vc £ • 
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Comparison between both energy scales gives a scaling estimate for the critical force density, 

Fc ~i (I) "? exp (^© 2 £)- 



(21) 



B. Dynamic approach 



In the dynamic approach [ |12|jl3| ] the critical force F c is determined directly without any reference to 
the pinning radius R c . To begin with, the system is driven by a large external force density F ex t S> F c . In 
this regime, the elastic manifold does not noticeably feel the pinning potential and F ext is only opposed 
by the friction force F f r i Ct which we describe by the phenomenological expression F f r i Ct = —770 v, with 
the viscosity rj a material constant. The flow velocity v of the manifold is determined by the steady- flow 
condition 

Ffrict + Fext = => v = v Q := — F ext ■ 

If F ext is decreased, the system begins to sensibly interact with the randomly distributed impurities via 
the pinning force F p i n . This becomes noticeable through fluctuations in the velocity field, leading to a 
reduction of the average velocity, v — v — <5v with v Q • 5v > 0. The average of F pin can be interpreted 
as an effective friction force with a velocity-dependent viscosity 5rj(v), 

(F pin )=-5r)(v)v. (22) 

As before, the steady-flow velocity is given by the relation 

F frict + (Fpin) + Fext = => r,(v) v := (tj + Sr](v)) V = F ext . (23) 

The applied force for which the relative fluctuations of the viscosity are of order unity then furnishes a 
useful criterion for the critical force density: 

Srj(F ex , 



= 1. (24) 

Alternatively, a criterion can be formulated making use of the reduction in the flow velocity, (F P i n ) = 
—r] S'v, and the condition 5v(v c ) = v c . 

We now calculate the friction coefficient Sr) due to pinning within a perturbative approach up to first 
order in 5 P . We start from Eq. (p2|), 

v a <fy = -(i^„(x,vi + u(x,*))), 

and expand the pinning force into a Taylor series, keeping the first two terms, 

v a Sr, ~ -(F«„(x, vt)) - (F p t(x, vt) ^(x, <)) . 

The first term vanishes; in the second one, we use the integral equation ( |l0| ) generalized to the time- 
dependent case and expanded to lowest order, 

W Q (x,i)= J d 4 yd S G(x-y,t- S )i^ m (y,v S ). 

Remembering Eq. (^3|) and the definition F pin = —E pin , one obtains 



v a 5r, 



J dtG{0,t)Kf\wt), (25) 



6 



where we made use of the time-translation symmetry. Inserting the explicit expression for the correlator 
Kf^P and for the time-dependent Green's function (|l^), we find: 



Srj(v) iV + 2/7r\§ . (°° dt ( v 2 t 2 \ ^ /r „ 



The integration produces a log-divergence for t — > 0. A lower cutoff is given by the time t% the elastic 
manifold needs to recover from a distortion on the smallest length scale £. Balancing elastic and dynamic 
terms in the Green's function G(x, i), Eq. (|l2|), on the scale £, we obtain 

(26) 

On the other hand, the exponential factor in the integrand, which originates from the correlator , 
provides an upper cutoff, 

^ ^ , (27) 

v 

describing the time scale for the interaction between the system and the pinning centers at velocity v. 
With these two cutoffs, the ratio between Sij and r/ becomes 



SyjFext) „ N + 2 
rj IQtt 2 



(§ 
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0(61). (28) 



This expression is valid in the non-linear regime, i.e., for forces F ext < 4v / 2C/£ (so that < t v ), where 
the effective friction coefficient rj(F ext ) = F ext /v(F ext ) deviates from rj . The critical force F c is found 
from the result (p8|), using the criterion Srj(F c )/r] ~ 1 (for weak disorder S p 1, we expect F c to be 
much smaller than C/£, so that the constant in Eq. ( |28| ) can be neglected). We obtain 

F c ~ -V- exp -— — - - , (29) 



Z \ N 

which is indeed exponentially small. As in the previous method, both the constant of proportionality and 
the numerical factor in the exponential function can only be ascertained up to a number of order unity. 



IV. RENORMALIZATION GROUP ANALYSIS 



As we have seen in the previous section, the leading term in the fluctuations of the displacement field 
u(x) shows a logarithmic behavior, (u 2 ) oc lni? 2 , cf. Eq. (|T9|). The logarithmic dependence weights all 
scales equally such that all length scales are relevant in the final result. This provides the motivation to 
apply the renormalization group (RG) . Our guiding line for this analysis will be the work of Efetov and 
Larkin fl4j| dealing with the pinning problem of charge- density waves. 



A. Construction of the RG transformation 

The following real-space renormalization procedure was suggested by Khmel'nitskii and Larkin |l5| ] for 
the problem of friction between two rough surfaces. In a recent publication B, Balents, Bouchaud, and 
Mezard gave a precise formulation of this method within the framework of the momentum-shell RG and 
used it to investigate the random manifold regime (large scales R > R c ) of the ((4-e)+AT)-model; here, 
we are interested in the perturbative regime (small scales R < R c ) with the aim to determine R c . 

Assume that we have renormalized the displacement up to a scale R\ > £; let us denote the corre- 
sponding field by Um(x). The associated free-energy functional reads 

-CVti£j(x) ■ Vuft)(x) +.Epi„ 1 (i)(x,U( 1 )(x)) . 
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Next, we go over to a larger scale i? 2 > Ri and separate um(x) into a far and a near-field contribution, 
U(!)(x) = U( 2 )(x) + w(x), with 

«f 2 )W = - / d 4 y^(x-y)B^ (1) (y,u (1) (y)), 

|x-y|>fl 2 

w a (x) = - [ d 4 j/G^(x-y)£;J„ (1) (y,u (1) (y)), 
where f2 := {x, y\Ri < |x — y| < i? 2 }- The free energy can then be written as 



(30) 



-CVu^ 2) (x) • Vuf 2) (x) + £ pi „ i(2) (x,u (2) (x)) 



with the renormalized pinning energy density, 

£ pi „,( 2 )(x,u (2) (x)) = £ pm , (1) (x,u (1) (x)) + -C [Vio a (x) • V M a (x) + 2 V«f 2) (x) ■ Vio°(x) 

The mixed term Vu?L • Vw™ is zero when integrated over and can be omitted. Integrating the quadratic 
term by parts, one finds 



^,m,(2)(x,u (2) (x)) = £p iri!(1) (x,u (1) (x)) - -w Q (x) CV 2 w a (x) 



(31) 



In the limit of zero temperature, out of all the possible near-field contributions w(x) only those minimizing 
the energy are relevant (in a statistical sense), 



dffpm,(2)(x,U(2)(*)) 

dw a 



= 0. 



Combining this condition with Eq. (j3l|), we can relate w(x) to the pinning force E™ in 

-CV 2 #(x)=-^ (1) (x,u^(x)). 



(32) 



(It is interesting to note that by differentiating Eq. ( pl| ) with respect to u%\> one shows that for the 
pinning force 



E pin,(2) ( X > U (2) ( x )) = E %in,{l) ( X > U (l) ( X )) + ° 



( dw 
V du 



(2) 



holds, which tells us that the force correlator K a,3 (u = 0) will not change under the RG.) 

Eq. (^p, together with Eqs. ( ^0| ) and (S3), is the starting point for our further considerations. It will 
allow us to derive a functional renormalization group (FRG) equation for the pinning energy correlator 
Kr(u), the subscript R denoting the scale of renormalization. We will then solve this equation by 
expanding the relevant correlators around u = 0. It turns out that the fourth derivative, K^f^ S {0), as 
well as higher-order even derivatives diverge at a finite scale which we identify with the collective pinning 
radius R c - Indeed, we know that on length scales smaller than R c (perturbative regime), the system is 
in a unique state and the pinning energy as well as the related correlators are analytic functions in the 
field u. At R c the manifold starts to probe different energy 'valleys' and our calculation, which rests on 
the analyticity of the pinning energy, breaks down, reflected by the singular behavior of the curvature of 
K(u). In what follows we will take the divergence of K%^ 5 (0) as an unambiguous definition of R c which 
allows one to go beyond the simple estimate (u 2 (R c )) ~ £ 2 . 



B. One-loop FRG equation 

We first construct a formal expansion in w of E pin (2), Eq. (|3~l|), which we will then use to derive a 
perturbation series in S p of (E pin ^ 2 )E pin ^}. Thereby, we will relate the energy correlators at renormal- 
ization scales i?2 and Kr 2 and K^, respectively; this will directly lead to the FRG equation for 
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Kr(u). To simplify the notation, we will omit the subscript 'pin' and use the abbreviations 



Ki. 



=- E 'pOT,(2)(Xj,U(2)(Xj)) , E\ 

: -^pin,(2)( x ji U (2)( X J')) i 
--Kr 2 (u (2) (Xj) - U' (2) (Xj)) , 



./a/3- 



a/3- 



: =^(l)( X J> U (2)( X j))> 
' :=£ 5n,(l)( X J» U (2)( X j)). 

:=^f> 2 ( Xi )-u 2 (x,)). 



Let us furthermore define the operator 

G?f(---) : / d^G^fc-xj) (--.). 
We begin by expanding the first term on the rhs. of Eq. ([n]) with respect to w(x) into a power series: 

oc ^ 

E (1) (xo, u (1) (x )) = (x , u (2) (x ) + w(x )) = |j S o 1 "' a ' (*o) ■ • • w Q! (x ) - 



i=0 

Next, we replace everywhere w by its implicit definition, Eq. (|30|). The force function E^ under the 
integral is in turn written as a Taylor expansion as above. We proceed iteratively, thereby eliminating w 
from our equation. Grouping the different terms with increasing number of Green's functions, we arrive 
at 

-25 (1) (xo, u (1) (x )) = -E + Eg E\ Eg Gft E\» G% E% - ±Eg» E% E° + 



+-^0 ^01 "12 ^2 ^23 ^3 + 21^0 "01 ^1 "12 ^2 ^13 ^3 



2 1 

1 7^M^ r~ i P jV T? V P (~<P a T? a r^ K ^ T?X 1 rpPPK f~iP>v t?v /~<p<J t?<t s-ik\ rpX 

2! 3 ^ 02 3 



1 

3f 



In a diagrammatic language this can be written as 



fj, v 



fi v p a 




2- 



(jt 




n 

1 



cr-2 

2 




fj, v p a k X 



A bullet with subscript 'j' and Z legs a,\, . . . , ai stands for the expression (— E" 1 '" 01 )/^., whereas a solid 
line with the indices '/i' and V joining two bullets H' and l j' represents G^" , i.e., the Green's function 
G M!y (xi — Xj), plus an integration over Xj £ f2. The factor 2 in front of the penultimate term is the 
symmetry factor which accounts for the fact that the expansion gives the two equivalent diagrams 



In a similar way, we expand the second term on the rhs. of Eq. (|31|), making use of Eq. ( |32| ) 
-i W «(x ) CV 2 w a (x Q ) = -i W "(x ) B ( a 13 (x 0) u (2) (x ) + w(x )) = . . . . 



Adding the series for i?(i)(xo, U(i)(xq)) and — iw a (xo) CV 2 w> a (xo) we finally obtain 




Let us now turn to the energy-energy correlation function, (EqE' ,), where xy :— x'y We write both 
energy functions in an expansion as done above. Using the linearity of the averaging procedure, we 
obtain a series of n-point correlators of the form (E^ ' ai ^ r " ■ ■ ■ E^ ) a "' 3 " ' ). Assuming Gaussian disorder, 
we can make use of Wick's theorem: n-point functions with odd n vanish and those with even n can be 
decomposed into a sum of products of 2-point correlators. For instance, application of Wick's theorem 
transforms the second-order term 

G^G^, (ESEiEjE'f) 

into the expression 

(E%Ej)(E»E°) + (EZE°)(E»EJ) + (EZE1){EJE'£) , 
diagrammatically represented by 



°01 ^O'l' 







1 





•- 



1 
-• 







1 

-m 



o' 



i' 



0' 



1' 



0' 



1' 



Dotted lines stand for 2-point correlation functions; the number of these lines indicates the order in 6 P 
of the corresponding diagram. 

The number of diagrams is considerably reduced by the special properties of the correlations. First, 
the correlator /sTr^u), being an even function in u, has all its odd derivatives vanishing at u = 0, 
K^" ai (Q) = (I = odd). Therefore, each term containing an odd correlator connecting two points T 
and l j' on the same 'tree' (i.e., already linked by a solid line) vanishes, 



■13, 



,^-^ 1 -/'n (u(2)(x . ) _ U(2)(x . )) 



= (-l) n K a R \ 



ai—aiPi—Pn(n\ _ 



(0) = (l + n = odd) . 



Furthermore, terms with a correlation between two neighboring points T and 'j' joined by a solid line, 
(Ef" Ej'") cx S 4 (xi — xj), give no contribution, because the delta function has no overlap with the 
domain of integration Q. An example is given by the last diagram above. The remaining non-vanishing 
diagrams up to second order in S p (zero and one-loop diagrams) are 
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• f • *r — 1* * — 1* * 

(E E' Qt )= \ + I \ \ + \ V +2-1 . (34) 



The factor 2 in front of the last term counts the number of equivalent contributions that appear when 
multiplying the two power series for the energies. Note that so far no special assumption about G M!y has 
been made. If the matrix of Green's functions is symmetric, G^ v = G"^, the second and third diagrams 
in Eq. (|34| ) are equal. 

The energy correlator K R2 (u(xa) — u'(x )) =: K 00 at scale i?2 is obtained by integrating Eq. ( p4| ) over 



K 00 = / d 4 x' (E E' ,). 

By way of illustration we discuss the calculation for the second diagram in Eq. (Q) , 

J d+x'o [D2] = -Kj» j d 4 Xl G^(x - xi ) G" ff (xo - x x ) Ktf . 

With R\ < |xo — xi| < i?2 and U( 2 ) smooth over we can replace by and extract it from 
under the integral. Inserting the explicit expression for the Green's function, Eq. (pd|), and going over to 
spherical coordinates, one easily carries out the remaining integral: 

S 3 [ R * dR 2 (R 2 

cTxi G(x — xij' = — 7— 

n 



where S3 = 2ir 2 is the surface of the unit sphere in 4 dimensions and / := S3/ (2(2ir) 4 C 2 ). The other 
diagrams are calculated in a similar way. Collecting all the terms, one finds 



K R2 (u)=K Rl (u)+I Q^(u)^(u)-^(u)^(0)) In 
which yields the functional RG equation for the pinning energy correlator, 



(36) 



= 1 (i ^p (u) ^p (u) „ ^p (u) AT(0) j ) (37) 

the initial condition being given by the function (u) , Eq. . This result was first obtained by Fisher 
for the (4+l)-model by means of the replica formalism. 



Eq. ( |37| ) is equivalent to a system of coupled differential equations for the expansion coefficients 
K^- ^O) (recall that all odd correlators K^'" aan+1 (0) = 0): 

d\nR 2 2 1Kr WKr (Uj ' 
dlnR 2 

^InR^ = 1 { K n 0tlP ^ K R W (°) + ^'""(O) + ^""(0) ^ w (0)) , (38) 

The first two coefficients are readily given by 



K R {0)=Ke{Q)(l--N8 p Jn tl 



I R 2 



(0) = ^" (0) = ^(0)1^, 
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where S a ^S a ^ = N and Eqs. (Oh and (Q) have been used. Eq. (38) can be simplified, if one assumes that 
the tensorial structure of K' R ^(0) is conserved under a RG transformation. In analogy with Eq. dl5|), 
we define Tr by K^f jS (0) =: TnA al3lS . Replacing this expression into Eq. (38), setting 7 = a and 8 = /?, 
and summing over all indices, one finds after some algebra 



dT 



R 



d\nR 2 



= (N + 8)ir 2 R =:A 1 T 2 R , (39) 



with the initial condition = 7u£ 4 (f ) 2 ■ This equation is easily integrated, 

E« = - (40) 

r e 1 - Ai in(i? 2 /e 2 ) ' 

the solution being valid for £ < i? < i? c . diverges at the collective pinning radius 

At that scale -ftTf^u) stops being an analytic function of u and for larger scales the system is no longer 
in a unique ground state. One may show by induction that all the higher derivatives, K^ 1 " Q2 ™ (0), n = 
3, 4, . . . , - if nonzero - exhibit a singularity at the same scale R c . The one-loop RG result for R c , Eq. (|4l|), 
confirms the lowest-order perturbation calculation, Eq. (po|), except for the argument of the exponential, 
which is smaller (for N < 8) by a factor 2JV/(Ar + 8). Note that because the renormalization is arbitrarily 
started at the length scale £, the present RG method cannot fix the constant of proportionality in Eq. (|4l|) 
just as the lowest-order perturbation calculation in section III. A. However, with this approach there is 
no unknown constant as in the condition (u 2 (R c )) ~ £ 2 ; the numerical factor inside the exponential 
function depends only on Re, i.e., the precise choice for the form factor, Eq. (^). Another advantage of 
this definition is that it provides a possibility for computing higher-order corrections to R c as will be 
done in the next paragraph. 

C. Two- loop FRG equation 

The number of diagrams inflates with increasing order of S p . There were only 3 non- vanishing second- 
order (one-loop) terms, there are already 20 non-vanishing third-order (two-loop) diagrams. With the 
assumption of symmetric Green's functions, = G vlL , only 14 topologically different diagrams are left. 
Two different integral expressions occur; they are of the form 

h= dxi G(x - xi) 2 / dx 2 G(xi - x 2 ) 2 , 

h = I dxi G(x - xi) 2 / dx 2 G(x - x 2 ) G(xi - x 2 ) , 
Jq[ Jn' 2 

where the integration domains are such that |xo — xi|, |xo — x 2 |, and |xi — x 2 | are between Ri and i? 2 . 
Accordingly, the two-loop diagrams can be regrouped into 2 classes; they are shown in Figs. 1 and 2. The 
factor in front of the diagrams enumerates the occurrence of this term in the series expansion. 
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Fig. 1: Two-loop diagrams in the class I2 



Fig. 2: Two- loop diagrams in the class I3. 

Recalling that a Z-leg bullet carries a factor we calculate these diagrams as demonstrated in the 
previous paragraph. The third-order contributions to Kr 2 , Eq. (|36|), are found to be 



K 



2 - l °°p(,,\ - 

R2 



(u) 



'3 ■ 



(42) 



Here, we have already taken into account the diagonal structure of the Green's matrix, — G. In 
section III.B we have seen that the relevant quantity for determining R c is Tr, which is essentially the 
fourth derivative of K(u). We therefore differentiate Eqs. (36) and j42| ) twice with respect to u a and 
twice with respect to vP ', set u = and use the definition of T R , K^ a0/3 (O) = T R A aat3/3 = T R N(N + 2). 
After summing over all indices, we eventually obtain 
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T Rl + T Rl (N + 8)h + T 2 Ri (N 2 + 6N + 20) J 2 + T 3 Ri (20N + 88) J; 



3 • 



Let us compute the two integrals I 2 and I3. I 2 is found at once, cf. Eq. ( p5|) 

h 

The calculation of I3 is less trivial. It can be written in the form 

h 



(I 1 ) 2 = I 2 In 2 ('|f 



(43) 



(44) 



{2^C 2 J ni x 



x 



d 4 y 



y 2 (x + y) 2 



where the integration domains f2i and 2 are such that |x|, |y|, and |x + y| are taken between R\ and 
R 2 . This integral is well known from </> 4 -field theory p6j. In the limit R\ — ► 0, R 2 — ► 00, it is 
logarithmically divergent, the divergent contributions being of the form ln(i? 2 /i?i) and ln 2 (i? 2 /i?i). To 
extract these terms, we take R 2 ^> so that the condition Ri < |x + y| < R 2 can be omitted with good 
approximation. We introduce spherical coordinates as usual, R := |x|,r := |y|, and choose the 2/4-axis 
parallel to x: 



47r5 3 
(2tt) 8 C 2 



^ dR?_ 
2~R2 



dr 2 



"1 



sin 2 dd 



R 2 



2Rr costf 



R 2 >i?i 



The dr 2 -integration is split into two parts, R\ < r 2 < R 2 and R 2 < r 2 < i? 2 ; for each part the dis- 
integration is carried out with help of the method of residua by transforming to coordinates z :— exp(ii9). 
The remaining integrals pose no problem and one obtains 



In I 12 



^ 2 

2 



R{ 



1 



m 



const. 



The change i?i — ► R 2 need not be infinitesimal. We are free to choose R 2 ^> Ri in such a way that 
ln(i? 2 /-Ri) 3> 1, as long as 5 p \n(R 2 /Ri) C 1, so that T R2 — T Rl remains small. For weak pinning, 
5 P <C 1, this condition is fulfilled and the last two terms in the above expression can be neglected, 



In 



i? 2 



Inserting the expressions for I 2 and is, Eqs. (|4j 
corrections up to two-loop approximation: 



1 + AlT Ri In 



A\T 2 Ri 



Im 2 
2 

and 
hi 2 



i?2 X 



i? 2 



(45) 



j) into Eq. (fl3|), we obtain the logarithmic 



R2 



A 2 T Rl In 



R 2 



(46) 



with Ai = (N + 8)1 and A 2 = (20iV + 88)/ 2 . The fact that the two-loop term cx In 2 is the exact square of 
the one-loop contribution (AiT Rl In) provides a check for the correctness of the one-loop RG equation, 
Eq. (39). Indeed, from Eq. (HO) one deduces 



r,B 2 



1 + AiT Rl In 



R 2 



A\Y 2 Ri In- 



i? 2 



0$) 



this shows that the In -contribution in Eq. 



is already taken into account in the one-loop equation 



and need therefore not be considered explicitly. The two-loop RG equation for Y R then reads: 

dT R 



A 2 T\. 



Similarly, the two- loop FRG equation for K R (u) is inferred from Eqs. ( pq ) and 
ln 2 -terms in I 2 and I3: 



(47) 

), by omitting all the 
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dK R (u) 
dhxR 2 



= I 



V K%>(u) - Kg(0) K£»(u) K£'(u 



(48) 



The integration of the differential equation (fl7|), 



in ir 



dlnR' 2 = 



In J 2 



(II 



yields the solution 



if) 



1 



r. 



1 



Ax r 2 + a 2 rs ' 

1 A 2 



dT 



An 



AxT 2 A 2 r A\{Y + A 1 /A 2 ) 
Ax 



Ax T R Ax T ? 



a 2 r R 



Ax 

a 2 r € 



To extract i? c , we take the limit Tr — > oo and use the assumption of weak pinning, so that A 2 T^/Ax oc 

i? 2 \ 1 A 9 



5 P < 1 



In 



i' 2 



In 



^2 £| 



Solving with respect to R c and inserting the explicit expressions for Ax, A 2l J, and Tj, leads to the final 
result for the collective pinning radius, 



RcZiZS™ exp 



ax 



(49) 



where 



ax 



N + 8 \ ir 



and a 2 — 



2{5N + 22) 
(N + 8) 2 



Comparison with the one-loop result, Eq. (|4l]), shows that the two-loop calculation produces an algebraic 
correction with exponent a 2 (the constant ax in the exponential function remains unchanged). Finally, 
scaling arguments provide us with an improved result for the critical force density, 



c ~ £ U 



£ P 



exp 



2ax 

8v 



D. Dynamic approach and RG 

In this paragraph we apply the above results to the dynamic approach introduced in Sec. III.B. We 
investigate the behavior of the friction coefficient rj(v) = rj +Sr](v) of Eq. ( p3[ ) under the renormalization 
group and arrive at an estimate for the critical force density F c in an alternative way. 

In order to pave the way for our real-space RG, we rewrite Eq. (|25|) through the Fourier representation 



of the Green's function, cf. Eq. (12) 



v 6r]= dtK^'^wt) I ——J- exp 

Jo 4 J (2tt) 4 V \ V 

where we have consistently replaced r] by r\. Since disorder is weak, the critical velocity v c is expected 
to be small. Hence, we expand the correlator K^^(vt) into a power series around wt — and retain the 
first non- vanishing term, proportional to T^: 

/»r „X F £ f , ( Ck2t 



15 



6r) = (N + 2)T i r 1 I J 



The time and angular integrations are easily carried out, 

k 2 ■ 

This form is appropriate to make contact to our results. Let us derive a RG equation for the viscosity: 
Assume we have renormalized rj up to a scale R\ > £. Going over to the larger scale R2, the viscosity is 
increased by 

VR2 ~ VRi = (N + 2)IT Rl r) Rl -jy, 

Jl/Rl « 



= (N + 2)IT Rl VRl ln^ 

As a result, the desired one-loop RG equation for the viscosity reads 

J^ = (N + 2) IT Rm . (50) 

The corresponding formula for the (4+l)-model has recently been derived by Leschhorn et al. 
using essentially the same method. Eq. (50) is solved by separation of variables, 

m dn f Vn d\nR 2 

^ = (N + 2)I dTT^^. (51) 

Use of the one-loop renormalization flow of T R , Eqs. ( j39|) and ([iol), yields 

i+ ^-?-(^)'-(* r « ta ir- <52 » 

with a = (A^ + 2)/(A^ + 8). 

In section III.B we defined the critical force F c as the applied force F ext for which 5-q(F ext ) ~ r) , 
reflecting the onset of nonlinearity in the velocity-force characteristics as F ex t is lowered. Here, we 
employ an alternative criterion, which we believe is better adapted to the problem. If a (finite) critical 
force F c > exists, then lim F t ^F + v (^ext) = 0. As a consequence, 

lim 7] (F ext ) = lim 



F cxt ^F+ F ext ->F+ v(F ext ) 

and we can define F c via the singularity of the effective viscosity rj = i] + 5r). Obviously, rj R diverges for 
R — > i? c , the critical force F c is thus related to the collective pinning radius R c . From scaling we know 
that 

C f 2 

F ext --^. (53) 



i r 2 

Recalling Eq. ©, wc thus find 

16tt 2 / 2 




iV + 8 \tt J S p 



(54) 



which essentially agrees with the lowest-order perturbation result Eq. ([29|), the argument of the expo- 
nential function being smaller by a factor a = (N + 2)/(N + 8). As usual, the number in front of the 
exponential remains unspecified. With the help of the scaling relation (|53"|), the viscosity, see Eq. (p52]), 
can be written as a function of F ex t, 



^t) = f ln F e .. 



16tt 2 \2J " \£F, 



(55) 



providing us with an explicit form of the velocity- force characteristics v = F ext / rj(F ext ) versus F ext . 
Whereas the lowest-order perturbation result, Eq. (pq), is only valid for applied forces F c <C F ext < 
C/£, this expression holds down to much smaller forces and thus better describes the non-linear regime. 
Expanding Eq. ( |55| ) for small 8 P and comparing the resulting series with Eq. one can check that 

both expressions coincide (for F ext €C/() up to first order in 5 P as they ought to. 
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V. SUMMARY AND CONCLUSION 



Using a RG approach, we have calculated the collective pinning radius R c as well as the critical force 
density F c for a (4+7V)-dimensional elastic manifold subject to point defects. We have confined ourselves 
to the case of zero temperature and have assumed weak and short-range disorder of Gaussian type. 

By means of a real-space renormalization procedure a functional RG equation for the pinning energy 
correlator Kr(u) was found in one and two-loop approximation. It has been shown that the fourth 
derivative of that correlator taken at zero u, K^f 1 (0) = r#A Q ' 37 ' 5 , possesses a singularity at some 
finite scale of renormalization which we have identified with the collective pinning radius R c . To one-loop 
correction we have derived the expression R c = a £ exp(otx / 6 P ) , a x = (2/n) N / 2 8tt 2 /(N+8), Eq. (|lj), that 
confirms the result obtained by a simple lowest-order perturbation calculation, Eq. (|2C|), the argument of 
the exponential function differing by a factor 2N/(N + 8). In two-loop approximation we have found an 
additional algebraic factor oc S^ 2 with an exponent a?2 = 2(5iV + 22)/{N + 8) 2 , see Eq. (^9|). Contrary 
to the perturbation calculation, the RG approach succeeds in fixing both ot\ and a,i. The constant of 
proportionality «o, however, remains undetermined in either of the two methods. 

Next, we have derived a RG equation for the viscosity r\ that appears when driving the system by an 
external force, cf. Eq. (pp[). The previously obtained result for Tr has allowed us to solve this equation, 
see Eq. (|5|). In comparison with the lowest-order perturbation result, Eq. (p8|), the RG solution provides 
us with an improved description of the velocity-force characteristics in the non-linear regime where the 
disorder potential considerably influences the vortex motion. We have defined the critical force F c via 
the singularity of the viscosity r\ and indicated how F c and R c are both related. As was the case for R c , 
the one-loop expression we have found for F c agrees with the result obtained by means of dimensional 
estimates and the dynamical approach, Eqs. (^lj) and (|29| ) respectively, except that the number in the 
exponential function is now fixed. 

The improved result for R c provides a firm starting point for the scaling analysis in the random manifold 
regime (large scales R > R c ). In a forthcoming publication pfj, we apply the same method to the physical 
case of a three-dimensional system of vortices in type-II superconductors with dispersive elastic moduli 
and gain quantitative refinements of the weak collective pinning theory. 
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